Generation of Dicke States with Phonon-Mediated Multi-level Stimulated Raman 

Adiabatic Passage 
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We generate half-excited symmetric Dicke states of two and four ions. We use multi-level stimu- 
lated Raman adiabatic passage (STIRAP) whose intermediate states are phonon Fock states. This 
process corresponds to the spin squeezing operation and half-excited Dicke states are generated dur- 
ing multi-level STIRAP. This method does not require local access for each ion or the preparation of 
phonon Fock states. Furthermore, it is robust since it is an adiabatic process. We evaluate the Dicke 
state using a witness operator and determine the upper and lower bounds of the fidelity without 
using full quantum tomography. 

PACS numbers: 



Genuine multipartite entanglement has been exten- 
sively investigated in quantum information science [l| . It 
is very difficult to generate genuine multipartite entangle- 
ment because the entangled states of many qubits decay 
with the square of the qubit number To date, the 
Greenberger-Horne-Zeilinger (GHZ) state of 14 qubits 
has been generated [lj using a M0lmer-S0rensen gate 
0; however, the half-excited symmetric Dicke states 
that form another important group of entangled states 
have been generated by photonic qubits with only four 
and six photons 0, 0] ■ No robust and efficient method 
has been demonstrated for generating Dicke states using 
more than two trapped ions[5|, [6j • 

Dicke states are generally the simultaneous eigenstates 
of the total angular momentum operator and its projec- 
tion along one direction. One subgroup of Dicke states is 
symmetric under permutation of qubits: 
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where Pk is the permutation operator and „C m = 
n\/[m\{N — m)!]. In this letter, we refer to these states 
as symmetric Dicke states. Here, we focus on the half- 
excited symmetric Dicke states (| D^,*)), which are dif- 
ficult to generate by addressing local qubits with an ar- 
ray of laser pulses [7[ orpost-selection by photodetection 
with atomic ensemble |8[. These multi-partite entangled 
states are not only important for quantum information 
science, they are also useful for precise measurements 
such as Heisenberg-limited measurements because 
these states are related to the squeezed vacuum state in 
quantum optics through the two-excitation correlation 
13]. Moreover, half-excited symmetric Dicke states can 
form decoherence-free entangled states with a dressing 
field 0. 

In this letter, we use phonon-meditated multi-level 
stimulated Raman adiabatic passage (STIRAP) 0,0 
to generate half-excited symmetric Dicke states of two 
and four ions. We evaluate the four-ion Dicke state with a 
witness operator and the fidelity without using full quan- 



tum tomography. 

We consider the case where ions are cooled to the 
ground state and are irradiated by two color laser beams 
that are detuned from the blue and red sideband tran- 
sitions for the M0lmer S0rensen interaction 0). In the 
interaction picture, the Hamiltonian is expressed by 
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where J± (= J] f a± , where cf± is the Pauli operator 
for the ith ion) is the global spin flip operator and a(a T ) 
is the annihilation (creation) operator for the center of 
motion of the ion chain, and S is the detuning of the blue 
and red sideband transitions. We begin with the initial 
state |44 • • • 4)« I 0) p , where | Q) p is the phonon vacuum 
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FIG. 1: (a) Energy-level diagram of phonon-mediated multi- 
level STIRAP. Vertical ladders indicate spin states and hori- 
zontal arrays represent phonon Fock states, (b) Evolution of 
spin noise for four ions. Solid, dashed, and dot-dashed curves 
indicate AJ^, AJy , and Ajf , respectively. 
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state. This initial state is obviously the dark state of this 
Hamiltonian when fif, = and we adiabatically change 
the amplitudes of the two color beams to transfer the 
population to the other dark states (|ft • • • t)s I 0) p , 
where Q r = 0). Here, we calculate the dark states and 
show that the half-excited symmetric Dicke states are the 
dark states of the Hamiltonian. 

We represent the Hamiltonian by the eigenstates of 
the global spin operator J 2 and J z . Due to the symme- 
try of the Hamiltonian, we can consider the symmetric 
subspace of the whole N qubit Hilbert space, whose ba- 
sis is | D™, >) a | n) p , where N is the number of ions and 
| n) p is the phonon Fock states. Fig. 1(a) shows the 
interaction for the lower phonon states (| 0) p , | l) p and 
1 2) p ). We assume that the red and blue sideband transi- 
tions have the same magnitudes for one-photon detuning 
5 and that the "vertical" Raman transitions [Fig. 1(a)] 
are on resonance relative to two-photon detuning. On the 
other hand, there are other Raman transitions, which are 
"horizontal" Raman transitions [Fig. 1(a)]; these transi- 
tions change the mean phonon number. For the case in 
which vertical Raman transitions are on resonance, hor- 
izontal Raman transitions have two-photon detuning as 
25. When 25 3> r]£l r ,r]fli,, the horizontal Raman transi- 
tions can be ignored and the interaction can be expressed 
in Fig. 1(a) corresponding to the multi-level Raman tran- 
sitions [6]. The Raman couplings are not independent 
and they have different strengths due to the cooperation 
effects of symmetric Dicke states. The coupling strengths 
of the cooperation effects are expressed by 
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and the Hamiltonian is expressed as follows in rotation 
frames 
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In this letter, we consider only even numbers of ions. 
This allows us to calculate the dark state | ipd(^r, ^b)) 
Hill] as 

N/2 
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where A is a normalizing constant. The dark state exists 
only for the motional ground state. 

Since each dark states (except for the initial and final 
states) is a superposition of states which are connected 
due to the quadratic raising and lowering operators of 
J+J+ and J_ J_, the dark states correspond to the vac- 
uum squeezed state in quantum optics IH 1(|. To see 
these characteristics, we calculate the spin noise for the 
four- ion case. Fig. 1 (b) shows the dependence of the vari- 
ance of the global spin operator (A J x , AJ y , A J z ) for dark 
states on the parameter 6 (fi(, = 1 — cos 6, fl r = 1 + cos 9) . 
For < < § or j < 9 < n, the spin noises are squeezed 
or anti-squeezed and when 9 — §(^6 = fi r )> spin noise in 
one direction goes to zero and the half-excited symmetric 
Dicke state in the x direction (| D 2 ^) = R y (ir/2) [D 2 ^)) 

is generated, where R v (0) is the global rotation operator 
about the y axis. Even for the general case of 2n ions, 
the half-excited symmetric Dicke state (iDJnfx))) ^ s S en_ 
erated when = Q r because the dark state [^(fi, f2)) 
is the eigenstate of the angular momentum operator in 
the x direction (J x = (J + + J_)/2) 

J x = o. 

This population transfer using the dark states of phonon- 
mediated multi-level Raman transitions is robust against 
certain fluctuations because the quantum state is always 
in the dark state of the Hamiltonian, as in the STIRAP 
method 
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We demonstrate phonon-mediated STIRAP with two 
and four calcium ions trapped in a linear Paul trap. De- 
tails of the experimental system used are described in 
our previous studies 0, [2(| . The center-of-mass mode 
of the axial direction is cooled by sideband cooling to the 
ground state (ncou ~ 0.02), while the other motional 
modes are cooled only by Doppler cooling. The qubit 
is composed of a ground state S 1( / 2 (™ = — 1/2) and a 
metastable state D 5 / 2 (ti' = —5/2) of a calcium ion and 
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FIG. 2: (a) Pulse sequences for phonon-mediated STIRAP. 
We measure the evolution of spin and spin noise by truncat- 
ing the STIRAP pulse after time r. (b) Time evolution of 
population transfer from the I-J444-) state to the |tttt) state. 
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we use S1/2 — D 5 / 2 transitions to excite sideband transi- 
tions. We first initialize all the qubits by optical pumping 
and irradiate them only with the red sideband laser pulse 
so that the initial state becomes the dark state. To oper- 
ate phonon-mediated STIRAP, we change the amplitude 
of the sideband pulses, as shown in Fig. 2(a). The peak 
Rabi frequency of transitions is r/fl = 2tt x 14kHz and 
the pulse length is t = 340/xs. The adiabatic condition 
is satisfied as t]Qt ~ 5 > 1 Fig. 2(b) shows the 

measured spin projection ( J z ) during population trans- 
fer for the four-ion case where a high-efficiency transfer is 
demonstrated. To see the characteristics of spin squeeze, 
the time evolution of the variance of the global spin oper- 
ator is measured. It is shown in Fig. 3 for four ions. This 
shows that spin noises are squeezed or anti-squeezed. 

When the two sideband pulses have the same ampli- 
tude, the half-excited symmetric Dicke state given by 
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is generated. We evaluate these half-excited symmetric 
Dicke states for the two- and four-ion cases. 

For the two-ion case, because the Dicke state is the 
same as one of the Bell states, we can estimate the fi- 
delity only with a global operation and by measuring a 
parity oscillation. When Q r = we irradiate the ions 
with a 7r/2 pulse whose phase is <f> and measure the ex- 
pectation value of a parity operator. Fig. 4(a) shows 
the parity oscillation; its amplitude is estimated to be 
A p = 0.95 ± 0.05 by fitting. Without the tt/2 pulse, 
the diagonal elements are measured to be {p-i,po,Pi} = 
{0.516, 0.033, 0.451} ± {0.011,0.004,0.011}. Using these 
values, we can obtain the fidelity as F = — 
0.96 ± 0.03, which is higher than 0.5 pH). 
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FIG. 3: (upper) Visualized spin noise in XY plane, which 
evolves with phonon-mediated STIRAP. The solid (dashed) 
curves shows J z < (J z > 0). (lower) Time evolution of spin 
noise. Circles, white squares, and black squares indicate AJ^, 
AJy and AJ Z , respectively. J x (J y and J z ) is (anti-)squeezed. 



To evaluate the four-ion Dicke state, we use (i) witness 
and (2) fidelity methods. For the four-ion Dicke states, 
the witness operator is given by [22[ 



Ji 



where i and j are any two orthogonal directions. For the 
Dicke state in the x direction (| D^^)), we can evaluate 

its entanglement with the witness W yz . After phonon- 
mediated STIRAP for four ions, we irradiate the ions 
with a 7r/2 pulse whose phase is <f> and measure the ex- 
pectation values of the square of the spin J| (Fig. 4(b)). 

This maximum peak (2.64 ± 0.05) corresponds to (Jy). 
On the other hand, ( ) is measured as 2.82 ±0.05 with- 
out an analysis pulse. The witness is then calculated as 

(W yz ) = 5.46 ± 0.07 > 5.23, 

which is greater than 3-a of the threshold (5.23) of gen- 
uine four-partite entanglement 22 1. 

We also evaluate the Dicke state with the fidelity, 
which is important to demonstrate that the generated 
state is closer to the Dicke state than the other genuine 
four-partite entangled states. For the Dicke state, the 
fidelity cannot be exactly determined without full quan- 
tum tomography }23j ] , w hich requires local access and 
many measurements [24|. Here, we evaluate the lower 
and upper bounds of the fidelity, rather than the fidelity, 
by using only global access and a few measurements. De- 
tails of this method are described in the Appendix. For 
the four-ion case, the fidelity(F) of the half-excited sym- 
metric Dicke state along the x direction is bounded such 
that 
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where p\ is the population of the angular 
momentum state whose spin component along 
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FIG. 4: (a) Parity oscillation of the Dicke state of two ions, 
(b) Squared spin of the half-excited Dicke state of four ions. 
The dashed line indicates shot noise for four ions. 
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the x direction is i. When O r = fit,, we 

(x) r (x) (x) (x) (x) 

measure p\ as {f>_2>P-i>f>o 1P1 iPa ; = 
{0.00, 0.03, 0.88, 0.03, 0.03} ± {0.00, 0.02, 0.03, 0.02, 0.02}. 
From these population and witness, we evaluate the 
fidelity as 

0.84 ±0.03 < F < 0.88 ±0.03 

Because the maximum fidelity between the GHZ and 
Dicke states is 3/4 (=| (GHZ | Dicke) | 2 , where |GHZ) = 
^{\mi)+ Itttt)), I Dicke) =| D 1(*)»' this lower bound 
(0.84 > 3/4) demonstrates that the generated state is a 
Dicke state. 

In summary, we demonstrate phonon-mediated STF 
RAP and generate the half-excited Dicke states of two 
and four ions. This method requires only the motional 
ground state and global access to whole ions. More- 
over, the adiabatic process, which corresponds to the spin 
squeezing operation, makes the method robust against 
certain disturbances and high-fidelity Dicke states can 
be generated. The fidelity of these Dicke states are mea- 
sured as F = 0.96 ± 0.03 for two ions and 0.84 ± 0.03 < 
F < 0.88 ±0.03 for four ions. Using the witness operator, 
we show that the generated four-ion state is a genuine 
four-qubit entangled state and, from the lower bound of 
the fidelity, we show that this state is the half-excited 
symmetric Dicke state. 

This work was supported by MEXT Kakenhi "Quan- 
tum Cybernetics" Project and the JSPS through its 
FIRST Program. One of the authors (N. A.) was sup- 
ported in part by the Japan Society for the Promotion 
of Science. 

Appendix 

We calculate the upper and lower bounds of the fidelity 
of the half-excited Dicke states (ID^J. We suppose that 
the witness {{W xy )) and the population of angular mo- 
mentum state along the z axis have been measured. 

W = (W X y) = (jj) + (jy) = (J 2 ) ~ (j|) 

p = {jP-miV-ju+U • • • iVOi • • • ,Pm} 

where jm is the maximum projection along the z axis, 
which is equal to the ion number divided by two. We con- 
sider the diagonal elements of the density matrix based 
on the angular momentum state along the z axis. 

Pj,j = (j,3z,Z I P I 3,3z,0, 

where | j,j z ,£) satisfies J 2 \j,j z ,£) = j(j + 1) \j,j z ,£) 
and J z | j,j x ,£) = j z | j,j z ,0- € (= l)2,--0 is an 
index that distinguishes between states with the same 
angular-momentum quantum numbers. States with j = 
]m comprise only these perfectly symmetric states. We 
then assume that the states | j — JM,jz,£ — 1) are 



symmetric Dicke states. We can calculate the population 
P such that pj z = J2j ■ £ Pj-jz-i- Note that Pj u ,o,i is the 
fidelity of the half-cxcitcd Dicke state. 
For the upper bound, we expand po as 

Po =^2Pj,o,£ ^ PiM.o.i = F 

j'.e 

because the diagonal elements of the density matrix are 
non- negative (pj,j s ,£ > 0). 

For the lower bound, we use the witness operator. As 
the witness operator is diagonalized by the angular mo- 
mentum state along the z axis, we have 

w = X)[7C7 + i)-jf]/»w.,e 

< j(J + l )PiV£ + Yl L?m(Jm + 1) - fflpj, ■ 

The first term on the right-hand side is again bounded 
as 

+ l)Wo,« 

< Jm(JM + 1)Pjm,Q,1 + 0'm - l)iM X! Pi>°>t 

< JmOm + l)Pj M ,o,i + Um - 1)j'm(po - Pmfi,i)- 

Using these inequalities, the fidelity has the following 
lower bound: 

W J M - 1 \- / ]M + l jz \ , F 

2^m 2 f-f V 2 2] M J 
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